Gravitational Waves from Rotating Neutron Stars 
and Evaluation of Fast Chirp Transform Techniques 


Tod E. Strohmayerf § 

f Laboratory for High Energy Astrohphysics, NASA/GSFC, Greenbelt, MD, 
20771 


Abstract. 

X-ray observations suggest that neutron stars in low mass X-ray binaries 
(LMXB) are rotating with frequencies from 300 - 600 Hz. These spin rates are 
significantly less than the break-up rates for essentially all realistic neutron star 
equations of state, suggesting that some process may limit the spin frequencies 
of accreting neutron stars to this range. If the accretion induced spin up 
torque is in equilibrium with gravitational radiation losses, these objects could be 
interesting sources of gravitational waves. I present a brief summary of current 
measurements of neutron star spins in LMXBs based on the observations of 
high-Q oscillations during thermonuclear bursts (so called “burst oscillations”). 
Further measurements of neutron star spins will be important in exploring the 
gravitational radiation hypothesis in more detail. To this end I also present a 
study of fast chirp transform (FCT) techniques as described by Jenet & Prince 
(2000) in the context of searching for the chirping signals observed during X-ray 
bursts. 


1. Introduction 

X-ray binaries are potentially among the most interesting sources of gravitational wave 
emission which current and future gravitational wave detectors will attempt to study. 
The high frequency signal produced during binary inspiral and ring down of black 
hole and neutron star binaries contains detailed information on the properties of the 
compact object as well as the structure of spacetime in its vicinity. These objects will 
be prime targets for ground based detectors such as LIGO which because of seismic 
noise are only sensitive in the high frequency range above ~ 100 Hz. 

Neutron stars are compelling targets of investigation because of the extreme 
physical conditions which exist in their interiors and immediate environs. For example, 
the signals produced by inspiral of a neutron star binary just before merger depend 
on the equation of state (EOS) at supranuclear density, a quantity which is still not 
precisely constrained by currently available astrophysics and nuclear physics data (see 
for example Heiselberg & Hjorth- Jensen 1999). Moreover, fundamental properties of 
the star, such as its mass, can be extracted if the gravitational wave signal can be 
measured. Thus gravitational wave astronomy can in principle provide new probes of 
fundamental physics as well as advancing neutron star astrophysics. 
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Table 1. Burst Oscillation Sources aild Properties 


Sources 

Frequency (Hz) 

Ai/ (kHz QPO, in Hz) 

References 1 

4U 1728-34 

363 

363 - 280 

1, 2, 3. 4, 5, 13, 14 

4U 1636-53 

290, 5SU 

251 

6, 7 

4U 1702-429 

m’ 

315 - 344 

4, 9 

KS 1731-260 

524 

260 

10, 11, 12 

Galactic Center 

589 

Unknown 

15 

Aql X-l 

549 

Unknown 

16, 17 

X1658-298 

567 

Unknown 

18 

4U 1916-053 

270 

290 - 348 

19, 20 

4U 1608-52 

619 

225 - 325 

8, 21 

SAX J1808-369 

401 

Unknown 

22,23 


’References: (1) Strohinayer et al. (1996); (2) Strohmayer, Zhang, k Swank (1997); (3) 
Mendez k van der Klis (1999); (4) Strohmayer k Markwardt (1999); (5) Strohmayer et 
al. (1998b); (6) Strohmayer et al. (1998a); (7) Miller (1999); (8) Mendez et al. (1998); 
(9) Markwardt, Strohmayer k Swank (1999) (10) Smith, Morgan, k Bradt (1997); (11) 
Wijnands k van der Klis (1998); (12) Muno et al. (2000); (13) van Straaten et al. (2000); 
(14) Franco (2000); (15) Strohmayer et al (1997); (16) Zhang et al. (1998); (17) Ford (1999); 
(18) Wijnands, Strohmayer k Franco (2000); (19) Boirin et al. (2000); (20) Galloway et al. 
(2000); (21) Chakrabarty (2000); (22) Heise (2000); (23) Ford (2000) 

2. Burst Oscillations 

Radio observations provided the first indications that some neutron stars are spinning 
with periods approaching 1.5 ms (Backer et al. 1982). These rapidly rotating neutron 
stars are observed as either isolated or binary radio pulsars. Binary evolution models 
indicate that neutron stars accreting mass from a companion can be spun up, or 
“recycled”, to millisecond periods (see for example Webbink, Rappaport & Savonije 
1983). This formation mechanism likely accounts for a substantial fraction of the 
observed population of millisecond radio pulsars, however, other formation scenarios 
have also been proposed (see van den Heuvel & Bitzaraki 1995; van Paradiis et al 
1997). 

In recent years direct evidence linking the formation of rapidly rotating neutron 
stars to accreting X-ray binaries has been provided by data from the Rossi X- 
ray Timing Explorer (RXTE). The first evidence came from the discovery of high 
amplitude, nearly coherent X-ray brightness oscillations (so called “burst oscillations”) 
during thermonuclear flashes from several neutron star binaries (see Strohmayer 2000 
for a recent review). These oscillations likely result from spin modulation of either one 
or a pair of antipodal “hot spots” generated as a result of the thermonuclear burning 
of matter accreted on the neutron star surface. Indisputable evidence that neutron 
stars in X-ray binaries can indeed be rotating rapidly then came with the discovery of 
the first accreting millisecond X-ray pulsar SAX J1808-369 (Wijnands & van der Klis 
1998; Chakrabarty & Morgan 1998), which is spinning at 401 Hz. 

The observed distribution of burst oscillation frequencies, including the 401 Hz 
pulsar, is very similar to the observed distribution of millisecond radio pulsars. Table 
1 lists properties of butst oscillations in the sources in which they have been detected. 
The RXTE observations suggest that the spin frequencies of neutron stars in accreting 
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Figure 1. Contour plot of the dynamic power spectrum of a burst from the 
LMXB 4U 1702-429 observed on April 1, 2001 (UTC) with the PCA onboard 
RXTE. Shown are contours of constant power spectral amplitude as a function of 
frequency and time. The burst profile is also plotted (right axis). The oscillation 
frequency increase from about 328 to 330 Hz in ~ 10 seconds. 


binaries span a relatively narrow range from 300 - 600 Hz. Moreover, these observed 
frequencies are significantly less than the maximum neutron star spin rates for almost 
all but the stiffest neutron star equations of state (Cook, Shapiro & Teukolsky 1994). 
In principle, accretion should drive the spin frequencies close to break-up, unless some 
other mechanism intervenes to remove the accreted angular momentum. 

3. Gravitational Wave Emission from Rotating Neutron Stars 

Bildsten (1998) recently proposed that the angular momentum gain from accretion 

could be offset by gravitational radiation losses if a misaligned quadrupole moment 
of order 10" 8 - 10~ 7 //vs could be sustained in the neutron star crust. Here I^s 
is the moment of inertia of the neutron star. In this scenario the strong spin 
frequency dependence of gravitational radiation losses sets the limit on the observed 
spin frequencies. Recent theoretical work has also shown that an r-mode instability in 
rotating neutron stars may also be important in limiting the spins of neutron stars via 
gravitational wave emission (see Bildsten 1998; Andersson, Kokkotas, & Stergioulas 
1999; Levin 1999; Ushomirsky, Bildsten & Cutler 2000; Andersson et al. 2000). 

The strong i/ 3 * * 6 frequency dependence of the energy radiated by gravitational waves 
means that rapidly rotating neutron stars with misaligned quadrupole moments might 
have observationally interesting gravitational wave amplitudes. The spin periods 
of neutron stars inferred from burst oscillations cluster rather tightly in the range 
from 300 - 600 Hz. White & Zhang (1997) suggested that the observed range of 
spin frequencies could be produced if these neutron stars were spinning in magnetic 
equilibrium. However, in order for the observed frequencies to be similar would require 
that mass accretion rate and the magnetic moment be correlated. It is not presently 
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known if such a correlation is to be expected based on theoretical grounds. The model 
proposed by Bildsten (1998) suggests that the spins of these neutron stars may be 
limited by the emission of gravitational waves (see also Wagoner 1984). The spin 
down torque due to gravitational wave emission is proportional to n 5 so that one 
would expect a critical spin frequency above which accretion torques are cancelled 
out by gravitational wave losses. By equating the characteristic accretion torque with 
the gravitational wave torque one can determine the average quadrupole moment 
required to maintain the critical frequency. For the mass accretion rates characteristic 
of LMXBs and a critical spin frequency of 300 Hz one obtains a quadrupole moment 
Q ~ 4.5 x 10 37 gm cm 2 , or about 5 x 10 -8 //vs (see Bildsten 1998). The question 
that remains is whether or not such a quadrupole moment can be routinely generated 
within a neutron star. This is presently an area of active research (see for example, 
Ushomirsky, Bildsten & Cutler 2000). 

Regardless of the mechanism, if the accretion torque is indeed balanced by 
gravitational wave losses then the amplitude of the gravitational radiation can be 
calculated (see Wagoner 1984; Bildsten 1998). The dimensionless strain h is in the 
range from h ~ 10 -27 - 10“ 26 . Although less than the estimated sensitivity for LIGO 
I, one can greatly improve the sensitivity by pulse folding if the rotational ephemeris 
of the neutron star is known (see Brady & Creighton 1999; Ushomirsky, Bildsten & 
Cutler 2000). Current estimates indicate that a narrow band configuration for LIGO-II 
will reach interesting search limits for these neutron stars, especially for the brightest 
of the LMXBs, for example, Sco X-l (Ushomirsky, Bildsten & Cutler 2000). This 
also provides strong motivation for additional deep X-ray timing searches in order to 
detect coherent pulsations in more LMXBs and to measure the pulse ephemerides so 
as to improve searches for gravitational wave emission. It also illustrates the strong 
synergism between X-ray and gravitational wave astronomy in the context of neutron 
star binaries. 


4. Searches for Chirping Signals: Burst Oscillations 

One way to search for additional neutron star spin periods is to try and make more 
sensitive searches for burst oscillations. Since the signals during many bursts are 
known to increase with time (ie. they are chirps, see Figure 1), it is likely that more 
general search techniques which allow for frequency changes will be better at detecting 
weak burst oscillation signals. Prince & Jenet (2000) have recently described a general 
method based on multidimensional FFTs to search for signals with varying frequency. 
They call this technique the Fast Chirp Transform (FCT). In an effort to develop 
more sensitive searches for burst oscillations as well as techniques to search for GW 
emission from rotating neutron stars I have begun evaluating the FCT on X-ray burst 
oscillation data which are known to contain chirping signals. Here we briefly describe 
the simplest 2-parameter (linear) chirp and the corresponding 2-d FCT. 

A signal with a linearly increasing frequency is of the form h a (t) = A(t ) sin(27r(i4-|- 
«>t 2 )), that is the phase advance is quadratic in time. The FCT is an approximation to 
the exact Fourier transform for such a signal; H koM = £\ o h jo exp(2Tri(k 0 (jo/N 0 ) + 

kiUo/No) 2 )). The second term in the exponential represents the quadratic phase 
advance produced by the changing frequency. The sum is over j 0 with 0 < j 0 < N 0 - 1. 
That is, the data, h jo , are a discretely sampled time series with N 0 bins. The FCT 
approximates this expression by breaking up the single sum into a double sum by 
defining a set of N\ subintervals of the data. Within each subinterval the phase 
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advance produced by the quadratic term is required to have the same value, that is, 
it increments uniformly as a “second index”, and the approximate expression can be 
evaluated as a 2-d FFT. This requires only O(N 0 Nilog 2 (N 0 Ni)) operations as opposed 
to O(Nq), and is therefore much more computationally efficient. There is a simple 
geometrical way to see how the 2-d FFT implementation of the FCT is constructed. 
One forms a two dimensional array, one dimension being N 0 in size (the number of 
samples in the original time series), and the other being of length AT, (the number 
of subintervals). The length of each of the N\ data segments packed along the N 0 
component is determined by the phase evolution function. Since the increase in phase 
is quadratic in time, subsequent data segments down the N i direction become shorter 
because they must represent the same total phase advance. In this way the boundaries 
of the subintervals are specified uniquely by the phase evolution function (see Prince 
& Jenet 2000). 

The range of i> searched is determined by N \ , the number of subintervals of the 
original time series. The maximum is searched is simply 0 max = N X /2T 2 , where 
T is the length of the original time series. Both positive and negative values are 
searched, analogously to the positive and negative frequencies of a standard FFT 
power spectrum, in this case, however, the FCT is not in general symmetric with 
regard to positive or negative i>. 

I have implemented the 2-d FCT in IDL using the formalism outlined in Jenet k. 
Prince (2000). I have evaluated the FCT by looking at burst oscillation data in bursts 
with detected oscillations from the LMXB 4U 1702-429 (see Markwardt, Strohmayer & 
Swank 1998). Figure 2 shows a comparison between a standard FFT power spectral 
analysis and one using the FCT for the burst shown in Figure 1. As can be seen 
in Figure 1, this burst has a strong oscillation at ~ 330 Hz. I used 8 seconds of 
data sampled at 4096 Hz (the interval between the dashed vertical lines in Figure 1). 
The data were first mixed with a f mix = 326 Hz modulation in order to bring the 
oscillation signal down to lower frequency, and therefore reduce N 0 . This is similar 
to a heterodyne procedure. Figure 2 (left) shows a plot of the peak power in the 
FCT as a function of is. One can see that the signal strength is greatly increased at 
is ~ 0.075 s -2 compared with is = 0. In Figure 2 (right) the power spectrum with 
the highest value of the FCT power is compared with that computed with is = 0 in 
Figure 2 (right). The FCT at positive is recovers a narrow, coherent peak whereas that 
computed with no frequency derivative is smeared out and weaker. This demonstrates 
the effectiveness of the FCT for detecting these kinds of signals. A detailed search for 
burst oscillations in new data will be presented elsewhere. 
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Figure 2. Fast chirp transform (FCT) analysis of a burst from 4U 1702-429. The 
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(dashed) and 0 — 0.078 (solid). 
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